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INTRODUCTION 


The  main  thrust  of  mathematical  analysis  of  deflagration  and 
detonation  waves  since  the  early  1970 's  has  been  through  the 
application  of  activation  energy  asymptotics  (ASA)  as  in 
Buckmaster  and  Ludford  (1982).  While  this  approach  is  of  great 
utility  in  many  problems  of  practical  interest,  there  are  also 
many  cases  where  the  activation  energy  cannot  be  taken  as  large, 
such  as  in  hydrogen  oxidation.  There  is  little  doubt  that 
hydrogen  will  be  an  important  fuel  in  the  coming  decades  since 
the  need  for  energy  sources  with  high  specific  energy 
(automobiles,  aircraft.  Space  Shuttle,  etc)  will  be  with  us  for 
the  foreseeable  future  and  concerns,  both  real  and  imagined,  over 
the  release  of  carbon  dioxide  into  the  environment  are  mounting. 

With  these  concerns  in  mind,  mathematical  methods  for  the 
study  of  combustion  phenomena  needed  to  be  developed  through  the 
application  of  large  heat  release  asymptotics  (LHRA) .  Previous 
work  in  this  area  for  deflagration  waves  is  scant  (Mikolaitis, 
1986a,  1986b,  1987a)  although  similar  adeas  have  been  used  in  the 
analysis  of  well  stirred  reactors  (Gray,  1973).  In  Mikolaitis 
(1986a),  the  structure  of  the  adiabatic,  plane  deflagration  wave 
with  one  step,  irreversible  kinetics  is  studied  under  the  limit 
of  large  heat  release.  In  Mikolaitis  (1986b),  a  similar  analysis 
is  undertaken  for  the  two-step  Zeldovich-Lifldn  model  that 
includes  a  chain  branching  mechanism  and  a  termination  step.  The 
motivation  of  this  work  was  to  show  that  the  usual  AEA  analyses 
of  this  problem  do  not  result  in  a  realistic  flame  structure  over 
the  entire  range  of  flame  temperatures  but  TjHRA  analyses  do. 

THE  ELEMENTS  OF  LARGE  HEAT  RELEASE  ASYMPTOTICS 

Rates  of  reaction  generally  follow  an  Arrhenius  rate  of  the 

form 
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rate  =  B  [F]3[0]‘=T^exp(-Ea/RT) 


where  [F]  and  [O]  are  the  concentrations  of  species  F  and  o,  T  is 
the  absolute  temperature,  R  is  the  universal  gas  constant,  B  is 
the  pre-exponential  constant,  a  and  b  are  the  stoichiometric 
integers  (these  need  not  be  integers  if  this  rate  expression  is  a 
model  for  a  more  complicated  reaction  mechanism) ,  n  is  the 
algebraic  temperature  dependence,  and  is  the  activation 
energy.  Eg/R  has  units  of  temperature  and  is  called  the 
activation  temperature. 

In  AEA  it  is  assumed  that  Ea/R  is  very  much  greater  than  the 
temperature  of  the  reacting  gases  everywhere.  A  consequence  of 
this  assumption  is  that  reaction  is  confined  to  a  thin  sheet  with 
frozen  reaction  on  one  side  and  equilibrium  on  the  other.  The 
diffusion  of  heat  and  mass,  on  the  other  hand,  take  place  on  a 
much  longer  length  scale. 

In  LHRA,  Ea/R  is  taken  to  be  of  the  same  order  of  magnitude 
as  the  flame  temperature  but  large  compared  to  the  reactant 
supply  temperature.  The  resulting  flame  structure  in  this  limit 
has  reaction  distributed  on  the  same  length  scale  as  the 
diffusion  of  mass  and  heat  and  so  the  reaction  zone  and  preheat 
zone  of  the  flame  are  of  similar  thicknesses. 

As  a  concrete  example,  let  us  analyze  the  structure  of  the 
downstream  far-field  of  an  adiabatic,  laminar  premixed  flame. 

The  adiabatic,  laminar  premixed  flame  problem  involves  the 
search  for  a  steady  solution  to  a  set  of  nonlinear  ordinary 
differential  equations  on  a  doubly  infinite  domain.  Cool,  fresh 
reactants  are  supplied  far  upstream  of  the  flame  and  hot,  burned 
gases  in  equilibrium  exist  far  behind  the  flame.  The  mass  flux 
of  material  normal  to  the  flame  is  a  parameter  in  the  set  of 
differential  equations  and  is  to  be  determined  as  part  of  the 
solution.  Determining  the  mass  flux  is  usually  the  goal  of  the 
analysis.  Of  course,  finite  differencing  for  the  numerical 
solution  of  the  governing  equations  subject  to  boundary 
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conditions  on  a  doubly  infinite  domain  is  not  a  practical  problem 
and  so  the  computational  domain  must  be  truncated.  , 

The  typical  computational  approaches  for  the  cold  boundary 
are  to  assume  a  zero  reaction  rate  between  x=-®  and  some  position 
x=a  and  integrate  the  resulting  equations  to  arrive  at  boundary 
conditions  at  x=a  or  simply  to  transfer  the  cold  remote 
conditions  directly  to  x=a.  Both  of  these  approaches  neglect  the 
"cold  boundary  difficulty".  Since  the  incoming  stream  of 
reactants  is  not  in  equilibrium  (otherwise  it  would  have  no 
chemical  enthalpy  that  could  be  turned  into  heat) ,  the  incoming 
stream  must  have  a  positive  rate  of  reaction  and  so  the 
conditions  at  infinity  are  not  a  solution  to  the  field  equations 
(the  reaction  term  is  positive  but  the  convective  and  diffusive 
terms  must  be  zero  and  so  cannot  balance  the  reaction  term) .  In 
practice  a  steady  solution  appears  to  exist  due  to  the  fact  that 
real  flames  only  propagate  finite  distances  and  the  rate  of 
reaction  in  reactant  stream  is  very  small.  This  allows  the  flame 
to  reach  the  combustible  mixture  long  before  any  appreciable 
reaction  takes  place  in  the  fresh  mixture  at  the  cool  reactant 
temperature.  One  of  the  triumphs  of  AEA  was  that  it  allowed  a' 
rational  analysis  of  the  steady  flame  problem  without  introducing 
in  an  ad  hoc  manner  any  extraneous  term  in  the  reaction  rate 
expression.  One  simply  realizes  that  the  "steady”  solution  is 
the  leading  order  solution  to  an  intrinsically  unsteady  problem. 
Here  we  have  generalized  this  idea  through  LHRA.  It  is  not 
necessary  that  the  activation  temperature  be  considered  as  large 
compared  to  the  flame  temperature  but  only  that  the  activation 
temperature  is  large  compared  to  the  cool  reactant  temperature. 
When  viewed  in  this  way,  AEA  is  seen  to  be  a  limiting  case  of  the 
mere  general  treatment  through  LHRA. 

For  the  hot  boundary  it  is  typical  to  specify  equilibrium 
conditions  at  some  position  x=a+L  where  L  is  "sufficiently" 
large.  The  sufficiency  is  usually  tested  by  increasing  L  and 
comparing  the  resulting  value  of  the  steady-state  mass  flux  with 
the  previously  calculated  value.  It  is  typical  to  consider  a  2% 
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change  as  small. 

Here  we  will  show  how  an  analysis  of  the  downstream  far- 
field  of  a  laminar  adiabatic  premixed  flame  can  be  accomplished 
through  a  series  of  simple  example  problems  that  are  of  interest 
in  relation  to  several  large  activation  energy  analyses  that  have 
appeared  in  the  recent  literature.  The  results  of  an  analysis 
such  as  this  can  be  used  to  apply  boundary  conditions  a  finite 
distance  behind  the  flame  consistent  with  the  far-field 
structure.  This  should  ameliorate  any  errors  introduced  through 
the  truncation  of  the  computational  domain.  It  may  also  make  it 
possible  to  use  a  smaller  computational  domain  in  practice  since 
the  boundary  conditions  derived  from  this  analysis  have  the 
correct  approach  to  equilibrium  built  into  them  and  so  this 
fairly  slow  process  need  not  be  reproduced  in  the  finite 
difference  scheme. 

There  are  no  conceptual  differences  between  the  examples  to 
be  investigated  here  and  examples  with  complete  chemical  kinetics 
and  sophisticated  multi -component  transport  models. 


A  FOUR-STEP  MODEL 

Changes  in  kinetic  mechanism  with  flame  temperature  have 
been  invoked  to  explain  certain  phenomena  peculiar  to  flame 
propagation  through  very  lean  hydrogen  mixtures  (Peters  and 
Smooke,  1985) .  The  Peters-Smooke  four-step  model  in  essence  is  a 
chain  branching  mechanism,  a  product  formation  reaction,  and  a 
reaction  that  consumes  chain  carriers  and  produces  an  inert 
intermediate  (chain  termination) .  At  sufficiently  low 
temperatures  the  chain  termination  mechanism  dominates  the  chain 
branching  mechanism  and  effectively  terminates  the  reaction  by 
eliminating  all  the  active  radicals.  Peters  and  Smooke  (1985) 
analyzed  the  response  of  a  premixed  flame  with  this  kinetic 
mechanism  to  small  aerodynamic  straining  (weak  flame  stretch) 
under  a  fast  recombination/large  activation  energy  limit.  Under 
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the  influence  of  weak  flame  stretch,  such  premixed  flames  with 
Lewis  numbers  less  than  one  (such  as  lean  hydrogen/air  mixtures) 
will  experience  a  rise  in  temperature  above  the  adiabatic  flame 
•temperature .  This  then  was  used  to  explain  how  hydrogen/air 
mixtures  nominally  below  the  lean  limit  can  support  cellular 
flames  and  why  the  troughs  of  such  flames,  where  the  flame 
experiences  positive  flame  stretch,  are  much  brighter  than  the 
crests  (which  may  be  extinguished)  where  the  flame  is  compressed 
and  hence  cooled. 

Although  this  explanation  is  very  appealing  (and  probably 
correct  for  the  most  part)  it  leaves  some  questions  unanswered. 
First  of  all,  an  0(1)  increase  in  flame  temperature  above  the 
adiabatic  value  is  required  to  move  the  flame  temperature  of  the 
very  lean  hydrogen-air  mixtures  above  the  transition  temperature. 
This  requires  0(1)  flame  stretch.  The  analysis  of  Peters  and 
Smooke  only  allowed  weak  flame  stretch.  In  addition,  HO2  is  not 
inert.  There  are  reaction  paths,  albeit  slow  ones,  through  which 
HO2  can  react  to  form  more  active  radicals  through  an  alternate 
branching  cycle  involving  the  formation  of  hydrogen  peroxide, 
H2O2.  Previous  analytical  work  in  detonation  wave  structure 
(Mikolaitis,  1987b)  has  shown  that  the  existence  of  this 
alternate  branching  cycle  has  an  0(1)  effect  on  the  length  of  the 
induction  zone  of  even  stoichiometric  hydrogen-air  detonation 
waves.  It  is  well  known  that  lean  premixed  hydrogen-air  flames 
show  appreciable  levels  of  hydrogen  peroxide.  It  is  also  well 
known  that  appreciable  changes  in  flame  speed  occur  if  HO2 
chemistry  is  ignored.  For  these  reasons  it  is  very  likely  that 
HO2  chemistry  plays  a  major  role  in  low  temperature  hydrogen-air 
premixed  flames.  Peters  and  Smooke 's  treatment  assumes  that  the 
slow  chemistry  of  this  species  is  infinitely  so. 

As  a  start  of  a  more  general  analysis  of  this  problem,  we 
will  first  introduce  a  four-step  reaction  scheme  that  has  both  a 
high  temperature  and  low  temperature  mechanism.  Secondly  we  will 
analyze  the  far  field  flame  structure  that  results  from  this 
mechanism  in  a  general  context  that  allows  0(1)  activation 
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energies  (as  opposed  to  restricting  the  validity  of  the  analysis 
to  the  large  activation  energy  limit) . 

REACTION  MODELING 

As  is  well  known  from  reaction  vessels  studies,  the  way  in 
which  oxygen  reacts  with  hydrogen  at  low  temperatures  is  very 
different  from  the  way  they  react  with  each  other  at  high 
temperatures.  The  source  of  this  change  in  kinetic  mechanism  is 
the  competition  between  the  reactions  H  +  O2  -►  0  +  OH  and  H  +  O2 
+  M  -*  HO2  +  M  where  M  is  a  third  body.  The  products  of  the  first 
reaction  readily  react  with  H2  to  form  additional  radicals  while 
HO2  reacts  with  H2  to  form  H2O2  (hydrogen  peroxide) .  The 
peroxide  intermediate  dissociates  into  two  hydroxyl  radicals 
through  a  relatively  slow,  large  activation  temperature  reaction. 
These  OH  radicals  then  quickly  react  with  H2  to  form  additional 
radicals. 

Here  we  will  study  a  simplified  kinetic  model  that  has  the 
following  features: 

1)  a  chain  branching  mechanism  that  dominates  at  high 
temperatures, 

2)  a  chain  branching  mechanism  that  dominates  at  low 
temperatures , 

and 

3)  a  product  formation/termination  reaction. 

This  can  be  most  easily  accomplished  through  the  four-step 
mechanism 

A  +  X  =  X  +  X  (1) 

X  +  X  »  P  +  P  (2) 

A  +  X  a  Y  +  Y  (3) 

X  +  Y  »  X  +  X.  (4) 

A  represents  the  major  reactants,  X  represents  the  highly 


reactive  radicals,  Y  represents  the  relatively  inert 
intermediates  and  P  represents  the  final  product.  Forward 
reaction  (1)  is  a  chain  branching  reaction,  (2)  is  the  product 
formation/termination  step,  (3)  competes  with  (1),  and  (3)  and 
(4)  taken  together  ((3)  +  2  x  (4))  is  another  chain  branching 
mechanism.  The  set  of  forward  reactions  (1-2)  forms  the 
Zeldovich-Liftd n  mechanism  that  has  been  extensively  studied  in 
recent  years  (Lifldn,  1971,  Seshadri  and  Peters,  1983,  Tam  and 
Ludford,  1984,  1985,  Joulin,  et  al,  1985,  Mikolaitis,  1986b). 

Numerical  subscripts  will  refer  to  the  reaction  numbers  1  to 
4.  Positive  subscripts  will  refer  to  the  forward  reactions  and 
negative  subscripts  will  refer  to  the  reverse  reactions.  All 
reaction  rates  will  be  assumed  to  be  proportional  to  the  mass 
fraction  of  each  reactant  and  have  a  general  temperature 
dependence.  For  example,  the  rate  of  reaction  for  the  reverse 
reaction  (1)  is  x2f_^(T).  Of  course,  Arrhenius  reaction  rates 
are  a  special  case  of  this  more  general  form. 


GOVERNING  EQUATIONS 

Under  the  assumption  of  calorif ically  perfect  gases  with 
constant  transport  properties  and  also  assuming  low  Mach  numbers, 
the  non-dimensional  governing  equations  describing  a  one¬ 
dimensional  reacting  front  with  the  aforementioned  four-step 
chemistry  are 

MT'  -  T"  +  qi(AXfi-x2f-i)  +  q2 (X2f2-p2f_2) 

+  q3(AXf3-Y2f_3)  +  q4 (XYf4-x2f .4 ) 

MA*  -  ^  A"  -  (AXfi-x2f_i)  -  (AXf3-Y2f_3) 

A 

MX'  «  X"  +  (AXfi-x2f_i)  -  2(x2f2-p2f_2) 

-  (AXf3-Y2f_3)  +  (XYf4-x2f_4)  (7) 


(5) 

(6) 
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(8) 


MY'  =  ^  Y"  +  2(AXf3-Y2f_3)  -  (XYf4-x2f_4) 

*  \  r 

I 

P  -  1  -  A  -  X  -  Y.  (9) 

M  is  the  mass  flux  normal  to  the  flame  (constant) ;  T  is  the 
temperature;  A,  X,  Y  and  P  are  the  mass  fractions  of  species  A, 
X,  Y  and  P  respectively;  Lx  and  Ly  are  the  appropriate  Lewis 

numbers;  and  the  q's  are  the  heat  releases.  Since  the  heat 
release  for  species  A  being  converted  into  species  B  must  be 
independent  of  reaction  path,  we  have  the  additional  relationship 

q4  =  (qi  -  q3)/2.  (lo) 

The  boundary  conditions  far  behind  the  flame  are  that  the  mixture 
is  in  equilibrium. 

THE  EQUILIBRIUM  STATE 

For  a  given  final  temperature,  T^,  where  the  subscript  b  - 
denotes  the  burned  state,  the  equilibrium  state  is  very  easily 
ascertained  in  closed  form.  If  we  define 


C  = 


-f4  +  y  f4  +  8(f_4  + 


ii!£2. 


1£i£z2 

^1^^3 


(11) 


and 

K  =  (f_^  +  a2  f_3)/(fi  +  f3)  (12) 

where  the  reaction  rate  functions  are  evaluated  at  the  burned 
temperature ,  T^ ,  then 

Ab  -  K  Xb,  (13) 
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Yb  =  c  Xb, 


(14) 


and  Xb  satisfies  the  quadratic 


Xb^ (-Kfi+f_i+2f2-2 (1+C+K) 2f_2+Kf3-c2f_3-Cf4+f_4) 

■h  4Xb(l-f-C+K)  f_2  -2f_2  =  0.  (15) 

The  state  of  the  burned  gases  is  linked  to  the  state  of  the 
fresh  reactants  through  the  total  enthalpy.  Multiplying  (6)  by 
(‘Il+q2/2)»  (7)  by  q2/2,  and  (8)  by  (qi-‘-qa2-q3)/2  and  adding  these 
results  to  (5)  yields  after  integration  and  the  application  of 
boundary  conditions 

Tf  +  (qi+q2/2)Af  +  q2Xf/2  +  (qi+q2-q3 ) ’jff/2  =  (16) 

Tb  (<3l'''q2/2)  Ab  +  q2^b/2  +  {qi+q2“q3)  ^b/2 

where  the  subscript  f  denotes  the  remote,  fresh  state. 
Specification  of  the  cool,  fresh  reactant  state  fixes  the 
constant  in  (16)  and  hence  substitution  of  (13-15)  into  (16) 
yields  a  single  equation  for  the  burned  gas  temperature. 


PHASE  SPACE  EQUATIONS  AND  PERTURBATIONS  TO  EQUILIBRIUM 

The  system  of  equations  (5-8)  are  autonomous  and  hence  the 
independent  variable  x  can  be  replaced  with  one  of  the  dependent 
variables.  In  most  cases  T  can  be  used  without  difficulty.  The 
system  is  thereby  reduced  in  order  by  one  and  the  domain  of  the 
adiabatic  plane  flame  problem  goes  from  being  doubly  infinite  (- 
«o<x<+«)  to  finite  (Tf<T<Tb  where  Tf  is  the  temperature  of  the 
cool,  fresh  reactants).  Letting  z  =  T',  the  governing  equations 
become 
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Mz  =  z  ^  +  qi(AXfi-x2f_i)  +  qs (X^f 2-p2f.2) 


+  q3(AXf3-Y2f_3)  +  04 ( XYf 4 -x2 f .4 )  (17) 


Mz 

dT 

1 

La 

'  dT^ 

-  fAXfi-x2f..j,  ) 

-  (AXf3-Y2f„j-) 

(18) 

MZ 

dX 

dT 

1 

Lx 

(z^) 

+  (AXfi-x2f_i) 

-  2(x2f2-p2f-2) 

- 

(AXf3- 

-Y2f.3 

)  +  (XYf4-x2f_4) 

(19) 

Mz 

dY  _ 
dT 

1 

Ly 

Z  ^ 
^  dT 

+  2 (AXf3-Y2f_3) 

-  (XYf4-x2f_4) 

(20) 

i 


along  with  (9) .  One  of  these  equations  may  be  replaced  by  an 
enthalpy  relation.  When  T=Tb,  equilibrium  holds  so  that 

A(Tb)  =  Ab,  X(Tb)  =  Xb,  Y(Tb)  =  Yb.  (21) 

and 

z(Tb)  =  0.  (22) 

The  price  that  one  pays  when  converting  a  problem  on  an 
infinite  domain  to  one  on  a  finite  domain  is  that  now  one  has  a 
boundary  value  problem  that  is  singular  at  both  endpoints.  Here 
we  will  be  numerically  integrating  the  governing  ecjuations  (17- 
20)  from  one  end  point  toward  the  other.  It  is  therefore 
necessary  to  analyze  the  local  structure  of  the  singularity  in 
order  to  specify  proper  initial  conditions  slightly  away  from  the 
singular  point. 

Perturbations  to  the  hot  singular  point  can  be  analyzed 
through  the  expansions 

z=«  f  ,  A=Ab+<^  »  X=Xb+<,x»  Y=Yb+<^f  T=Tb+«»’,  «<<1.  (23) 

and  the  boundary  conditions  (21-22)  become 

r(0)=  (0)=x{0)=^(0)=0.  (24) 
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The  non-linear  limit  equations  resulting  from  the  expansions  (23) 
are  easily  solved  subject  to  the  initial  conditions  (24) .  The 
system  admits  a  solution  of  the  form 

5  =ar  ,  ,  f-gr  (25) 


that  satisfy  (24)  if  a,  b,  c  and  g  satisfy  the  algebraic 
relationships 

Mg  =  g2  +  q^Ai  +  q2A2  +  q3A3  +  q4A4 
1  -i 

Mga  =  —  g'^a  -  A^  -  A3 
A 

Mgb  =  ^  g2b  +  2A3  -  A4 


Mgc  =  ^  g^c  + 


Ai  “  2A2  “  A3  +  A4 


where 

Ai  =  aXfafi  +  AfaCfi  +  AbXbfi  -  2XbCf_i  -  Xb^fii 
A2  =  2XbCf2  +  Xb^fi  -  (1-Ab-Xb-Yb:2fi2 
+  2(1-Ab-Xb-Yb) (a+b+c)f_2 
A3  =  AbCf3  +  Xbaf3  +  AbXbfi  -  2Ybbf_3  -  Yb2fi3 
A4  =  YbCf4  +  Xbbf4  +  XbYbf^  -  2XbCf_4  -  Xb2fl4. 


(26) 

(27) 

(28) 
(29) 


(30) 

(31) 

(32) 

(33) 


The  temperature  dependency  functions  (f^  functions)  and  their 
derivatives  are  to  be  evaluated  at  the  burned  gas  temperature  Tb- 
The  derivative  of  z  with  respect  to  T  is  identical  with 
T"/T'.  Since  T  approaches  a  constant,  Tb,  as  x  goes  to  infinity, 
T”/T'  is  negative  for  large  values  of  x.  This  implies  that  dz/dT 
must  be  negative  in  the  neighborhood  of  Tb  and  so  g  cannot  be 
positive. 

A  solution  to  this  system  of  algebraic  equations  can  be 
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found  with  g=0  and  a,  b  and  c  satisfying  a  set  of  three  linear 
algebraic  equations  whose  solution  in  closed  form  is  easily 
accomplished.  The  implication  of  g=0  in  the  neighborhood  of  T-T5 
is  that  the  conditions  of  the  burned  gas  far  behind  the  flame  is 
dominated  solely  by  chemical  reaction  with  convection  and 
conduction/diffusion  playing  no  role.  This  spatially  homogeneous 
solution  also  corresponds  to  the  final  approach  to  equilibrium  of 
a  quiescent  gas  in  an  adiabatic  enclosure.  Of  the  multiple 
solutions  to  the  set  (26-29),  this  is  not  the  appropriate  choice 
as  can  be  easily  determined. 

Letting  a=qi+q2/2,  ^=q2/2,  7=(qi+q2“q3)/2 ,  it  is  easily 
shown  that 

M(T  +  aA  +  +  7Y)  =  Z(1  +  +  C  (34) 

is  true  globally  where  C  is  a  constant  that  can  be  evaluated  at 
either  the  hot  or  cold  singular  point.  If  (34)  is  differentiated 
with  respect  to  T  and  evaluated  at  T=Tb  with  the  further 
restriction  that  dz/dT  =  0  at  that  point  (so  that  g=0) ,  then 

1  +  aa  +  fic  +  7b  =  0.  (35) 

This  constraint  is  inconsistent  with  the  solution  to  the  three 
linear  algebraic  equations  for  a,  b  and  c,  i.e.  (27-29)  with  g=0, 
unless  the  three  Lewis  numbers  are  all  exactly  one.  The  case 
where  all  the  Lewis  numbers  are  unity  will  be  discussed  later. 
Therefore,  although  there  is  a  solution  for  g=0,  that  solution 
violates  a  global  constraint  and  hence  is  of  no  interest  in  this 
problem.  We  shall  now  consider  only  solutions  with  g  strictly 
less  than  zero. 

Solutions  to  (26-29)  with  g<0  can  be  found  numerically.  The 
computational  difficulty  is  significantly  reduced  through  the 

derivation  of  the  auxiliary  equation 

_ M(1  +  aa  +  Be  +  Tfb) 

1+  (aa/LA)  +  (fiC/Lx)  +  {7b/Ly) 


g  = 


(36) 


to  be  used  in  place  of  one  of  the  original  equations.  Through 
the  use  of  this  equation  not  only  is  the  unwanted  root  with  g=0 
avoided,  the  permissible  solution  space  in  the  abc  space  is 
easily  located  between  the  planes  where  the  numerator  and  the 
denominator  of  (36)  are  zero.  Since  in  practice  the  Lewis 
numbers  L^,  Lx  and  Ly  can  be  quite  close  to  one,  it  is  clear  that 
the  permissible  solution  space  can  become  two  very  narrow  wedges. 
In  such  situations,  knowledge  of  these  bounds  on  the  solution 
space  is  critical  in  determining  a  suitable  "first  guess"  for 
numerical  root  finding  on  this  system  of  non-linear  algebraic 
equations. 


THE  SPECIAL  CASE  OF  UNIT  LEWIS  NUMBERS 

When  all  the  Lewis  numbers  are  identically  one,  there  is  a 
further  reduction  in  the  complexity  of  the  problem.  Equation 
(34)  can  be  integrated  once  again  to  give 

T  +  aA  +  +  7Y  =  constant  (37) 


so  that 


1  +  0 


(38) 


We  therefore  now  know  one  of  the  a,b,c's  explicitly  in  terms  of 
the  other  two.  Furthermore,  g  can  now  be  found  in  terms  of  two 
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of  the  a,b,c's  through  the  quadratic  solution  of  (26),  that  is 


g  =  (M  -  yM2-4(2Aiqi))/2,  (39) 

where  the  negative  root  is  taken  as  the  only  possible  way  to 
generate  a  negative  value  of  g. 


THE  REVERSIBLE  ZELDOVICH-LINAN  MECHANISM 

A  subset  of  the  reaction  mechanism  (1-4)  consisting  of  the 
forward  reactions  1  and  2  has  received  considerable  attention 
through  the  years  (Zeldovich,  1948,  LiA4n,  1971,  Seshadri  and 
Peters,  1983,  Tam  and  Ludford,  1984,  1985,  Joulin  et  al,  1985, 
Mikolaitis,  1986b) .  If  the  alternate  branching  mechanism  is  very 
slow  compared  to  the  main  chain  branching  mechanism  and  it  is 
found  that  there  is  negligible  amounts  of  Y  at  equilibrium,  then 
the  low  temperature  mechanism  (3-4)  can  be  ignored.  This  is  the 
typical  case  when  the  flame  temperature  is  high.  Here  we  will 
analyze  the  two-step  mechanism  (1-2)  where  each  reaction  is 
reversible  in  order  to  take  into  account  equilibrium  dissociation 
of  the  product.  This  is  not  an  element  of  the  earlier  analyses 
of  this  mechanism. 

The  equilibrium  state  as  a  function  of  burned  gas 
temperature  is  easily  determined  to  be 

Xb  —  /f^+yf  /f  '  ™  ^  *  ^“^b”^b  (40) 

— 1  1  2  —2 

where  the  temperature  dependency  functions,  f^,  are  evaluated  at 
the  burned  gas  temperature  T^}.  Expanding  A,  X  and  T  in  the 
neighborhood  of  the  equilibrium  point  through 

A  =  Ab+«J,  X»Xb+<x»  T  =  Tb+er,  Z  =  ef  (41) 

leads  to  a  set  of  nonlinear  ordinary  differential  equations 
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subject  to  the  initial  conditions 


(r-0)=x-(r=0)=f  (r=0)=0.  (42) 

This  set  of  nonlinear  differential  e  ions  admit  solutions  of 
the  form 


S=ar,  r=cr,  f:=gr  (43) 

if  a,  c  and  g  satisfy  the  algebraic  equations 

Mg  =  g2  +  q^Bi  +  q2B2  (44) 

Mga  =  —  g2a  -  Bi  (45) 

La 

Mgc  =  ^  g2c  +  Bi  -  2  B2  (46) 


where 

Bi  -  aXbfi  +  AfeCfi  +  AfaXbfi'  -  2XbCf_i  -  Xb^f-i'  (47) 

and 

B2  =  XbCf2  +  Xb2f2'  +  2(1-Ab-Xb) (a+c)f_2 

-(1-Ab-Xb)2f_2'.  (48) 

These  solutions  also  satisfy  the  initial  conditions  (42) .  As 
before,  there  is  a  solution  to  this  set  of  equations  with  g=0, 
but  again  this  solution  is  not  appropriate.  In  order  to  avoid 
this  unwanted  solution  if  numerical  root  finding  is  used  it  is 
advantageous  to  rearrange  (44-46)  to  arrive  at 


^  Md-t-aa+flc) 

^  l+aa/L^+^C/Lx 


(49) 


provided  that  not  both  and  Lx  are  unity. 

If  both  Lewis  numbers  are  unity  then  an  exact  solution  with 
physically  sensible  behavior  is  easily  achieved.  In  that  case  a 
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and  c  are  related  through 


1  H  aa  +  =  0.  (50) 

(45)  is  then  divided  by  a  and  subtracted  from  (44) .  Elimination 
of  c  from  the  resultant  yields  a  quadratic  in  a. 


CONCLUDING  REMARKS  ON  THE  EXAMPLE  PROBLEM 

The  approach  to  equilibrium  behind  the  adiabatic  premixed 
flame  is  a  very  slow  process  that  typically  takes  place  on  a  long 
length  scale  compared  to  the  primary  reaction  zone  thickness  and 
the  heat  conduction  -  major  species  diffusion  length  scale.  Here 
we  have  analyzed  the  structure  of  the  approach  to  equilibrium 
behind  a  premixed  flame  front  for  a  four  step  reaction  model . 

The  simultaneous  non-linear  ordinary  differential  equations  that 
describe  the  process  can  be  solved  exactly  up  to  a  set  of 
constant  coefficients.  The  required  coefficients  are  found 
through  a  set  of  simultaneous  algebraic  equations  that  can  be  ' 
solved  numerically. 

The  mathematical  algorithm  used  in  this  example  problem  can 
also  be  used  in  the  general  context  of  a  complete  chemical 
kinetic  mechanism  with  variable  transport  coefficients.  We  are 
currently  investigating  the  ozone  decomposition  flame  as  an 
example  problem.  Through  such  a  solution  the  boundary  conditions 
an  infinite  distance  behind  the  reacting  front  can  be  brought  to 
a  finite  distance  behind  the  reacting  front  without  a  substantial 
loss  of  accuracy.  It  is  also  possible  to  reduce  the 
comp\J.tational  domain  for  the  calculation  of  an  adiabatic  plane 
premixed  flame  structure  significantly  because  the  long  approach 
to  equilibrium  will  be  known  analytically  and  need  no  longer  be 
resolved  numerically. 
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A  NEW  SIMILARITY  SOLUTION  FOR  REACTING  GAS  FLOWS 


In  the  AEA  analysis  of  premixed  flames  the  analog  to  the 
stagnation  point  flow  similarity  solution  is  very  similar  to  the 
non-reacting  stagnation  point  flow  in  that  the  gas  velocities  at 
infinity  are  singular.  Since  this  is  the  case,  the  applicability 
of  the  analysis  is  restricted  to  the  boundary  layer  applications 
where  the  flame  is  restricted  to  be  close  to  the  stagnation  plane 
even  though  the  solution  is  a  solution  of  the  full  Navier-Stokes 
equations. 

Under  the  limit  of  LHRA,  this  is  no  longer  the  case. 
Solutions  for  reacting  opposed  flow  are  found  where  the  velocity 
at  infinity  is  finite  and  uniform  without  slip  as  in  boundary 
conditions  for  flow  exiting  a  porous  plug.  This  means  that  the 
realistic  boundary  conditions  are  met  for  the  full  problem  and 
restriction  to  a  boundary  layer  solution  is  not  necessary.  A 
copy  of  the  page  proofs  of  an  article  on  this  solution  is 
included  as  Appendix  A. 

Why  is  the  character  of  the  solution  different  in  the  limit 
of  large  heat  release?  The  essential  difference  between  the  two 
solutions  is  that  in  the  AEA  stagnation  point  flow  type  solution 
there  is  a  quadratic  variation  in  pressure  transverse  to  the 
flame  whereas  in  the  LHRA  solution  the  transverse  pressure 
gradient  is  required  to  be  zero  to  leading  order.  There  is  a 
mathematical  argument  to  show  that  this  must  be  true,  but  a  non¬ 
rigorous  physical  argument  can  be  made  that  is  more  revealing. 

In  the  study  of  gaseous  bubbles  rising  through  liquids, 
where  there  is  a  great  difference  between  the  density  of  the 
gaseous  phase  and  the  liquid  phase,  the  density  ratio  between  the 
liquid  and  gas  phases  can  be  approximated  by  infinity  (i.e.,  zero 
density  for  the  gas  phase  to  leading  order) .  Under  these 
circumstances  free  streamline  theory  can  be  applied  and  in  this 
theory  the  material  surface  separating  the  high  momentum  outer 
flow  of  the  liquid  phase  from  the  stagnant  (to  leading  order) 
wake  fluid  in  the  gaseous  phase  is  a  surface  of  constant 
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pressure.  In  the  LHRA  analysis  of  opposed  flow  we  have  a  very 
similar  situation  where  the  flame  plays  a  role  analogous  to  that 
played  by  the  material  surface  in  free  streamline  theory.  It 
should  therefore  not  be  surprising  that  LHRA  predicts  zero 
transverse  pressure  gradient  and  low  momentum  far  behind  the 
flame. 

The  LHRA  analysis  also  predicts  some  very  interesting 
behavior  that  cannot  be  predicted  through  AEA.  It  is  well  known 
that  when  premixed  flames  with  Lewis  numbers  less  than  one  are 
positively  stretch,  the  flame  temperature  rises.  If  adiabatic 
conditions  are  enforced  far  behind  the  flame,  the  AEA  analysis  of 
such  a  stretched  premixed  flame  would  never  predict  extinction 
and  the  flame  temperature  would  continuously  rise  with  increasing 
stretch.  If  conductive  heat  loss  to  the  gases  far  behind  the 
flame  is  allowed  by  specifying  a  low  temperature  for  the  remote, 
burned  products,  then  as  the  flame  stretch  is  increased  the  flame 
temperature  will  eventually  fall  (but  only  after  the  flame  has 
retreated  far  into  the  opposed  flow)  and  the  flame  structure  has 
become  quite  diffuse.  This  mathematical  behavior  has  been 
interpreted  as  extinction  even  though  mathematical  flame 
structures  still  exist.  In  LHRA,  a  new  extinction  behavior  is 
found  that  is  very  easy  to  understand  and  remarkable  in  that  it 
at  first  seems  counter-intuitive:  Flames  can  extinguish  by 
becoming  too  hot!  The  physics  of  this  mechanism  is  very  simple. 
In  Arrhenius  kinetics  the  rate  of  reaction  has  a  temperature 
dependence  that  goes  monotonically  from  0  to  a  positive  constant 
as  the  temperature  goes  from  zero  to  infinity  and  the  approach  to 
zero  is  through  exponentially  small  functions.  The  rate  of 
reaction  is  proportional  to  this  quantity  times  the  concentration 
of  the  reactants.  These  concentration  go  monotonically  from 
Infinity  to  zero  as  the  temperature  goes  from  0  to  infinity  and 
the  approach  to  Infinite  concentrations  is  through  algebraically 
large  functions.  This  means  that  the  rate  of  reaction  is  zero 
for  both  0  and  infinite  temperature  and  finite  and  positive  in 
between.  The  rate  of  reaction  therefore  has  a  maximum  at  some 


18 


well  defined  temperature.  The  LHRA  analysis  of  the  far-field  of 
stretched  premixed  flames  shows  that  a  necessary  condition  for 
flames  to  exist  is  that  the  flame  temperature  is  at  or  below  the 
temperature  that  gives  the  maximum  reaction  rate.  Since  positive 
flame  stretch  increases  the  flame  temperature  in  flames  with 
Lewis  numbers  less  than  one,  an  upper  bound  on  the  flame  stretch 
that  such  a  flame  can  endure  is  found. 

It  should  be  pointed  out  that  the  similarity  solution  given 
in  Appendix  A  is  a  much  closer  analog  to  experiment  than  the 
stagnation  point  flow  solutions  whenever  the  flame  is  away  from 
the  stagnation  plane. 


CONTINUING  WORK  IN  LHRA 

In  the  short  term,  work  is  continuing  in  LHRA  along  two 
fronts:  the  analysis  of  ozone  decomposition  flames  and  the 
generation  of  other  premixed  flame  structures  with  multi-step 
kinetics.  The  ozone  flame  is  a  particularly  attractive  example 
problem  because  it  has  been  studied  through  conventional 
nvmerical  approaches  (Heimerl  and  Coffee,  1980)  and  through 
conventional  AEA  methods  (Rogg  and  Wichman,  1985,  amoung  others) . 
The  results  of  this  analysis  (undertaken  as  part  of  a  Master's 
program  by  a  graduate  student  at  the  University  of  Florida, 

Department  of  Aerospace  Engineering,  Mechanics  and  Engineering 
Science)  should  be  more  accurate  than  both  the  numerical 
solutions  where  2%  convergence  was  deemed  adequate  and  the  AEA 
solutions  where  reaction  is  mathematically  restricted  to  a  much 
thinner  region  than  is  seen  in  practice. 

In  the  long  term,  it  is  hoped  that  this  research  will  result 
in  faster  computing  methods  for  the  prediction  of  flame 
structures  for  both  stretched  and  unstretched  premixed  flames 
through  the  generalization  to  complete  chemical  kineteics  of  the 
mathematical  results  generated  so  far.  It  is  envisioned  that 
this  process  can  be  completely  automated  so  that  the  user  need 
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simply  enter  the  kinetic  scheme  that  is  to  be  used  and  the 
program  will  automatically  construct  the  far-field  solution  and 
so  forth. 

Another  application  of  LKKA  is  in  the  study  of  diffusion 
flames.  This  will  require  a  different  methodology  because  the 
temperature  is  not  monotonic  through  the  flame  structure.  There 
are  other  variables  that  are,  however,  and  so  this  is  not  an 
insurmountable  problem. 
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On  the  Abrupt  Extinction  of  Premixed  Flames  with  Lewis 
Numbers  Less  Than  One 

DAVID  W  MIKOLAlTir,  i.)i‘nnrtn)riil  nt  Ai;rGsp,ici;  Etiairiahnnrj  Mi}chnnic:<.  and 
Eii(/in^nr,nif  Scinni:!::  tjmvnrsiiy  of  Florida.  Giiinesville.  Fl  32611 

(  < Max  22  /v.s’^  Iff  ittuti  fttfm  Jutu-  27.  /V-S.**) 

Abstract — Mere  uc  h;ivc  ;inaly/e(l  the  emcrniiig  equations  of  rcjcting  low  Xtach  number  (low  for  plane 
premised  llamc  slruciurcs  in  non-uiiiform  llotv  liclds  given  that  the  activation  lemperaiure  of  the  iniviure 
IS  Urge  compared  !o  the  lemperaiure  of  the  remote  reaclanis.  An  extinction  mechaniMYi  is  unciivercil  lor 
flames  With  Lc\m-.  numbcr^  leN>  than  one  that  cannot  he  found  through  the  u%e  t>f  activation  encre\ 
as>nip1(Mics.  a  limit  vvhere  the  activatutn  icm|Krature  of  Utc  mixture  i>  large  compared  to  the  burned  e.is 
icmpcraturc. 


INTRODUCTION 

A  central  problem  in  combustion  theory  is  to  understand  the  influence  of  non- 
uniform  flous  on  burning  rate.  A  particularly  simple  model  problem  that  has  been 
C-Xiensivcly  studied  through  the  use  of  activation  cncrg>  asymptotics  (.AEA)  is  the 
premixed  flame  in  an  opposed  flow  vvith  .A  —  B  irreversible  kinetics.  Analysis  of  this 
problem  have  been  performed  under  the  constant  density  approximation  (Buckmuste.  f 
1979.  and  Buckmastcr  and  Mikolaitis.  1982)  and  with  full  fluid  mechanical  coupling  ' 
under  a  boundary  layer  approximation  (Libby  and  Williams.  1982.  1983.  1984,  and 
Libby.  Linan.  and  Williams.  1983).  None  of  these  studies  have  revealed  an  abrupt 
extinction  mechanism  for  adiabatic  premixed  flames  with  Lewis  numbers  sufficiently 
less  than  one. 

Here  we  will  perform  an  analysis  of  the  premixed  flame  in  an  opposed  flow  without 
making  a  boundary  layer  approximation  or  assuming  that  the  activation  energy  is 
large.  All  we  will  assume  is  that  the  activation  temperature  is  large  compared  to  the 
temperature  of  the  cool  remote  reactants  and  that  the  burned  gas  temperature  is  on 
the  order  of  magnitude  of  the  activation  temperature.  In  this  way  the  concept  of  an 
adiabatic  flame  speed  is  retained  by  making  the  reaction  rate  of  the  cool  reactants 
exponentially  small  compared  to  the  reaction  rate  in  the  flame  itself.  Reaction  will  no 
longer  confined  to  a  thin  sheet,  however.  Instead,  reaction  will  be  distributed  on  a 
length  scale  comparable  to  the  length  scale  of  the  temperature  field.  This  has  proven 
to  be  a  more  realistic  flame  structure  for  relatively  hot  ozone  flames  (Mikolaitis.  1987) 
and  in  all  likelihood  will  apply  to  many  other  situations. 

This  is  also  in  accord  with  the  results  of  CotTcc,  Koilar.  and  Miller  ( 1984)  who  have 
studied  the  overall  reaction  concept  as  applied  to  50%  ozone-50%  diatomic  oxygen. 

40*/«  hydrogen-60y<i  air,  and  stoichiometric  mcihane-air  premixed  flames.  At  nomi¬ 
nal  laboratory  conditions  for  the  stale  of  the  cool  reactants  ( .'00  K  and  1  atmosphere), 
the  global  activation  temperatures  arc  3974  K.  3271,  and  1 5  680  K  with  burned  gas 
temperatures  of  1943  K,  2215  K,  and  2232  K.  respectively.  These  give  burned  gas 
temperature  to  activation  temperature  ratios  of  0.489.  0.677.  and  0.143  for  56"/« 
ozone-50%  air.  40%  hydrogcn-60%  air.  and  stoichiometric  mcihane-air  premixed 
flames.  Such  numbers  will  be  considered  as  0(\).  The  corresponding  ratios  of  cool 
reactant  temperatures  to  the  activation  temperature  arc  0.0755.  0.0917.  and  (1.0191. 

These  will  be  considered  as  small. 


<.1 


1 


/»  H  MIKWI  \/(ls 


U'o  will  iii'iidmicDMoniili/c  our  ijovcriimL’  cijii.iliou-.  in  the  n.imic  in.iiini’i  .is  in 
Mikol.iilis  ( l')S7)  wliori;  ;i  Minil.ir  problem  w.is  .inniv/eil  iiiulei  ilic  resirn  non  ol  iiiiii 
1  ess  IS  ninnbei  As  1  csvis  n umber  etVeels  b.ive  proven  ui  be  sery  uiiporl.iiil  in  llie  A I  A 
;in.!lv-.esol  llnmes  ssith  lUiw  non- uniformiiies.  ii  svouUI  seem  likely  ihnl  the  snnie  ssoniil 
.ipply  here  Siiii  e  wc  are  imercsied  m  low  speed  ilcdaeraJion  ss  uses,  a  loss  M.ieli  minitiei 
lipproximation  is  .ippropriatc  so  ibal  the  resiiicii!;  sleails  :s*,o  dinien .ii-i;.i;  eosernin.; 


<‘(m/);('v  -1-  i'((jj  )  (  r  - 

0 

(1) 

liijii  l  \i'\  +  flwVT.V'r  =  V- T  i-  Qij 

/■"  >■  exp  (IT) 

ill 

<'(/>// >’)  ('.V  -1-  ('(w>')/VV  =  liLV'y  — 

oT’' )'  exp  t  -  1/D 

(.■<) 

umi'm  f  v  ovru/Py  =  —Ppidx  +  P(V'ii  -t- 

(i'(ru/Px  +  i'v/rt)/r.x)/y\ 

(4) 

om'v  ('.V  +  <jvrv/t'y  =  —  rp/py  +  P(  V'v  -h 

U(f'//.rx  +  t'v/i'v)  I'  f)/?) 

(5) 

/>  =  1  r 

(f'l 

lor  eontinuiiy.  energy,  species  conservation,  v  and  y  momentum,  and  equation  of 
stale.  Wc  base  assumed  in  the  derivation  of  these  equations  caloriftcally  perfect  gases 
of  equal  specific  heats  and  constant  transport  properties.  The  parameters  L  and  P  arc 
the  Lewis  number  and  the  Prandtl  number  respectively.  T  is  the  temperature  non- 
dimensionalized  sviih  respect  to  the  activation  temperature  of  the  reaction,  u  and  v  arc 
the  -v  and  f  components  of  velocity. /» is  the  non-dimensionalizcd  variation  in  pressure 
about  the  ambient  combustion  pressure,  and  K  is  the  mass  fraction  of  reactants. 
Lengths  have  been  scaled  so  that  the  coctficient  of  the  Laplacian  of  T  is  unity  and 
velocities  have  been  scaled  so  that  the  coefficient  of  the  reaction  term  is  unity.  Q  is  the 
non-dimensional  heat  release.  The  details  of  the  non-dimensionalization  can  be  found 
in  Mikolaitis  ( I9S7(. 

Plane  flames  w  here  the  .v-axis  is  normal  to  the  flame  and  the  i  -axis  is  parallel  to  the 
name  (sec  Figure  I)  will  require  that 

T  =  T(.x).  Y  =  }'(.v).  «  =  i/(.v).  (7) 

The  equation  of  state  immediately  gives  g  =  q{x)  so  that  we  can  define 

Af  =  g{x)u(x)  =  Mix).  (8) 

Continuity  now  requires  that 

V  =  -vTdMldx.  (9) 

This  .V  momentum  equation  (4)  now  shows  that  cpjvx  is  a  function  of  .v  only.  This  then 
implies  that 


c'pjPxcy  =  0.  (10) 

Substituting  the  right  hand  side  of  (9)  for  v  and  A/  for  qu  in  the  r  momentum 
Fq,  (5)  gives 

»>/n-  =  -\(Pd-(TdMldx)ldx-  ~  SUiTilM'tixVdx  +  TitlMltl.x)').  (II) 
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Applying  the  condition  (10)  results  in 


*  P(i'(Til.\f!d\)ldx-  -  .UtHTdM;d\)dx  +  TUIM  d\)'  =  const  (12) 

£  so  in  cenerul  one  can  expect  quadratic  variations  in  pressure  in  the  r  direction 

r  (transverse  to  the  flame)  for  plane  flames  as  long  as  the  constant  in  ( 1 2)  is  non-/ero. 

Later  we  will  demonstrate  that  this  quadratic  variation  is  inconsistent  «ith  the  large 
heat  release  limit  so  that  the  constant  must  be  zero. 

The  energy  and  species  conservation  equations  (2-3)  become 

.UdTIdx  =  d-TIdx-  +  Qr-'  ycxp(-i,ri  (I3) 

MdYIdx  =  (!//.)</=  r/r/.\-  -  r-'  Te.xp  (- l/T).  (14) 

We  now  have  a  svstem  of  3  non-linear  ordinary  diflcrcntial  equations  to  solve  given 
by  (12-14). 


LARGE  HEAT  RELEASE  ASYMPTOTICS 

The  central  notion  of  the  approximation  scheme  that  we  will  be  using  it  that  the 
temperature  of  the  remote  upstream  reactants  is  small  compared  to  the  activation 
temperature  of  the  reaction.  Since  the  activation  temperature  is  also  the  temperature 
scale,  this  then  requires  that  the  remote  upstream  temperature  is  zero  to  leading  order. 
Since  the  remote  burned  state  must  be  a  singular  point  of  the  system  the  reaction  rate 
there  must  be  zero.  This  can  be  accomplished  by  taking  gases  far  downstream  as  fully 
burned.  In  order  that  the  "regular”  solution  is  found,  we  w  ill  specify  that  far  upstream 
of  the  flame  the  balance  in  the  energy  equation  is  between  conduction  and  convection 
with  reaction  being  much  smaller.  This  will  prevent  any  "fast"  deflagration  wave 
solutions  (see  Kassoy,  1985,  for  example,  for  a  discussion  of  fast  deflagration  waves, 
that  is  deflagration  waves  where  the  remote  upstream  balance  in  the  cold  reactant  gas 
stream  is  between  convection  and  reaction)  from  appearing.  As  a  result,  some  of  the 
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lv)tiiiil;ir\  cmuliliiiiis  ,iic 


/  (  -  » 


(1.  )■(  /  )  =  r  ,.)'(+/)  =  0.  iH'  (/a  -•  0  ;is  V 

(</’7  '</a')  (A/f/7'(/\)  —  I  ;i.s  A  —  -  /  . 


+  y  . 
(15) 


IJororc  ivodistuss  !hc  rest  (it'ihc  boiiiuiiiry  coiulilions  ncccss;iry  lo  produce :i  unique 
soliilion.  xvc  will  invcsligulc  how  the  boundary  condiiioiv.  used  so  far  intiuenee  the 
ehoieo  of  inieeraiion  eonsiani  in  (12).  As  a  —  —  the  mass  (lux  M  must  approach 
(</'  T/J\  )  (i/iT,  </a)  and  so  for  larjte  ncgatisc  v 

PiTiT’iDT  -  (T(T-iry)ir-ir)  +  tht'ity)'  ^  k.  (I6) 

where  ■  denotes  derivatives  with  respect  to  .v  and  K  is  the  constant  of  integration  in 
(12).  If  any  of  the  three  terms  approaches  a  non-zero  constant  as  .v  —  -  r..  then  alt 
the  temts  become  of  the  same  order  of  magnitude  and  T  grows  quarticully.  In  fact, 
an  exact  solution  to  (16)  is  given  by  T  =  AV/(27  —  6P).  From  this  we  can  get  the 
rest  of  the  unknowns  (//.  v.  />.  a)  for  an  exact  solution  to  a  viscous,  conducting 
incompressible  (low  problem. 

The  only  other  way  to  have  the  left  hand  side  approach  a  constant  is  if  two  of  the 
terms  become  very  large  but  their  sum  approaches  a  constant.  This  becomes  impos¬ 
sible  because  if  two  of  the  terms  become  of  the  same  order  of  magnitude,  then  it  is 
easily  determined  that  T must  attain  quartic  growth  as  before.  This  then  implies  that 
the  only  permissible  value  of  the  constant  is  zero  and  so  now  our  set  of  three 
governing  equations  is  (12-14)  with  the  constant  for  equation  (12)  being  0.  This 
condition  implies  that  there  is  no  transverse  pressure  gradient. 

The  physical  reason  why  there  cannot  be  a  non-zero  transverse  pressure  gradi¬ 
ent  is  that  the  density  ratio  of  the  hot  product  gas  stream  to  the  cold  reactant 
stream  becomes  infinite  in  the  large  heat  release  limit.  Therefore  the  interface  between 
the  hot  and  cold  gases  behaves  in  a  manner  similar  to  the  liquid-gas  interface 
of  a  rising  bubble  {i.e..  free  streamline  theory).  A  practical  combustion  problem 
that  may  have  this  character  is  a  premixed  llame  attached  to  a  bluff  body.  In  that 
configuration  the  flame  sheet  roughly  separates  cold,  high  momentum  Ruid  in 
the  free  stream  from  hot.  low  momentum  fluid  in  the  recirculation  zone.  Since 
the  density  difference  between  the  frecstream  and  the  wake  can  be  large,  say  a 
factor  of  7  or  8.  the  separation  streamline  might  be  profitably  considered  as  a  free 
streamline.  Since  the  wake  gasses  arc  so  slowly  moving  in  the  recirculation  zone 
of  an  attached  premixed  flame,  it  may  be  better  to  analyze  the  flame  structure 
through  the  limit  given  here  rather  than  through  a  boundary  layer  analysis  where 
the  speed  of  the  remote  reactants  is  taken  to  be  very  much  larger  than  the  local 
flame  speed. 


THE  REMOTE  BURNED  STATE 

An  analysis  of  appropriate  to  the  remote  burned  state  can  also  be  performed.  In  this 
case  we  will  transform  the  problem  into  a  phase  space  with  T  as  the  independent 
variable  through  the  definition 


;  =  r/r  </v 


(17) 
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III  .irriM;  .il  llic  m.1  hI 

/’:  i/(r</( /-■</'/(//  )  <//  )•<//  “  M:  (l{l z  list UH  )  ill  +  Iz  iil.M  tIT)'  =  (t  (ISj 
W  =  ilzitir  +  Q  T'  '  Y  c\p  {  -  \ir\:  (19) 

M  ,n  ill  -  O.UiKzilYnlTMiH  -  r  '  )i;xp(-l  /  )  r  (2'l) 

aiul  ihc  partial  >ct  oC  boundary  conditions  (IS)  become 

>■(0)  =  >•  ,.  YiT,)  =  (I.  r(0)  =  0.  r(r,)  =  (). 

(f/r/r/D/iV/ -*  I  as  r-O  (21) 


where  1]  is  the  final  temperature  of  the  burned  gases.  T,  is  to  be  determined  as  part 
of  the  solution. 

We  will  now  look  for  solutions  to  this  system  in  the  neighborhood  of  the  point 
r  =  T),  )■  =  0.  r  =  0.  M  =  M„  through  the  leading  order  c.xpansions 

T  =  T,  +  i:t.  =  r.:.  Y  =  t:Y.  U  =  M„  ~  .>»/.  r.  4  I  (22) 

appropriate  to  this  problem. 

The  system  governing  the  behavior  in  the  far  field  is 

ii;,(h  =  M„  -  Q  T;  '  .f  e,\p  (-IT,  l,:  (23) 

iH;ilylth)(h  =  L(iV/..r/f,r/r  +  T,"-' .1  e\p  (- IT,)  0  (24) 

il'ni  th'  =  —  3(</C/’</T)(r/’«»'</r')/C  —  (d'ldx)'  (dm  ciz)  —  (d'^'dr)(d/»idzYi 

4-  M„((dy  dr)(dm/dT)/i'  -f-  (d'midi')  C)  P  (25) 

with  initial  conditions 

;(0)  =  .1(0)  =  Mt(0)  =  0.  (26) 

The  general  solution  to  this  underdetermined  non-linear  system  is 

C  =  .4t  (27) 

.1-  =  Bx  (28) 

m  =  Cltj'  (29) 

where 

A  =  (LM„  -  ^(L'Ml  +  4Z.-.  ))  2  (30) 

B  =  M„(\  -  DA  Qy  -  L,Q  (31) 
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1  =  A/„  I/’  (32) 

■/  =  n  '  c-xp  ( -  1/7; )  (33) 

aiul  ('  is  arbitrary.  Casting  the  solution  h  ick  into  the  original  variables  we  then  get 
the  approximate  solution 

r  5:  .(Cr  -  /;)  (.34) 

K  s:  B(T  -  T,)  (35) 

M  M„  +  C{T,  -  T)'  (36) 

if  we  take 

6  =  r;’.  (37) 

A  number  of  points  can  be  addressed  with  this  solution  for  the  far  field  structure. 
If  the  solution  for  the  perturbation  to  the  mass  flux,  hi.  is  to  make  any  sense  we  must 
have  a  >  0  or  C  =  0.  In  the  case  where  C  =  0  we  can  take  6  as  an  arbitrarily  small 
scale  compared  to  e.  This  then  implies  that  the  variation  in  M  within  the  r;  neighbor¬ 
hood  of  r  =  r,  is  identically  zero  and  so  every  derivative  of  M  with  respect  to  T  must 
vanish  at  T  =  Tf.  We  therefore  know  that  .1/  must  be  identically  zero  in  some  finite 
neighborhood  of  r  =  7}or  else  it  has  an  essentially  singularity  at  T  =  T,.  Numerical 
calculations  on  the  full  Eqs.  (18-20)  indicate  the  former  and  furthermore  show  that 
M  is  identically  constant  from  T  =  0  to  T  =  7).  Of  course  the  boundary  conditions 
at  r  =  0  are  not  met  in  general,  but  for  a  specific  choice  of  Af„  they  arc.  That  choice 
is  the  adiabatic  burning  rate. 

Since  our  interest  is  in  the  effects  of  flow  non-uniformities,  and  therefore  varying 
M.  attention  must  be  placed  on  the  case  where  C  /  0.  When  this  happens  we  must 
conclude  that  is  negative  in  order  that  1  is  positive.  The  reason  for  this  is  that  as 
.V  -•  00  either  T'  is  positive  and  T’  is  negative  of  T’  is  negative  and  T’  is  positive  if 
T approaches  a  constant  value.  This  implies  that  dztdT (which  is  the  same  as  T'jT') 
must  be  negative  at  T  =  7).  We  therefore  know  that  .4  is  negative.  Of  course,  the 
Prandtl  number  P  is  positive  and  hence  if  i  is  to  be  positive  then  Mo  must  be  negative. 
This  result  is  not  too  surprising.  A  more  interesting  conclusion  can  be  reached  once 
the  condition  of  non-negative  mass  fractions  is  applied.  In  the  far  field  the  imposition 
of  non-negative  mass  fraction  requires  that  5^0.  This  will  be  true  for  any  negative 
value  of  M,  if  L  ^  I  but  will  not  be  true  in  general  if  7.  <  I .  When  7.  <  1  a  critical 
mass  flux  +  x  given  by 

Moc  =  -  L)  (38) 


can  be  uncovered.  A  vanishes  at  this  particular  value  of  My.  If  Mn  is  smaller  than  Moc. 
then  B  as  defined  by  (31)  becomes  positive  and  hence  the  mass  fraction  of  reactant 
becomes  negative  in  the  far  field.  Since  this  is  impossible,  we  must  have 

-V7.y/(l  -  M.,  <  0  (39) 


when  7,  <  I. 
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The  jneqiKilily  (3‘i)  can  be  rewriUen  in  ihe  (orm 

7;  '  exp  <  -  1/7;)  ^  (I  -  /.)  A/,i  L  0.  (40) 

!f  II  >  1.  then  (he  Icl'l  hand  ^idc  of  {40)  i\  a  i’.v!n  '.mnic  illi  .nc: easing  funclio:'.  of 
[,  ihai  increases  w  ithout  bound,  l  or  a  given  value  of  A/„  and  Lew  is  number  le.ss  ih.iii 
one.  a  minimum  possible  burned  g.i.s  lemperaiurc  is  predicted.  This  situation  is  of 
limited  interest  because  it  corresponds  to  a  reaction  rate  that  increases  without  bound 
with  temperature.  This  is  hardly  realistic  in  thai  lor  sulliciently  higli  temperatures  the 
fraction  of  collisions  that  have  the  requisite  amount  of  energy  approaches  one  while 
the  density  of  the  reactants  approaches  zero.  Wc  thcrelbrc  know  that  the  rate  of 
rwiction  must  go  to  zero  as  the  temperature  goes  to  infinity. 

In  the  ease  where  n  =  1 .  again  a  minimum  burned  gas  temperature  is  predicted  but 
the  left  hand  side  of  (40)  no  longer  grows  without  bound.  Since  the  left  hand  side 
of  (40)  is  bounded  between  0  and  !  for  all  positive  values  of  T,  .  A/,;  is  bounded 
between  0  and  7.;(1  —  L).  In  this  ease  a  minimum  value  of  \(„  is  predicted,  namely 
—  JLlil  —  L).  that  is  attained  at  an  infinite  flame  temperature.  At  any  finite  flame 
temperature  the  value  of  must  be  higher.  Again  this  case  is  of  limited  interest 
because  now  the  rate  of  reaction  goes  to  a  positive  constant  instead  of  zero  as  the 
temperature  goes  to  infinity. 

When  II  <  I  the  left  hand  side  of  (40)  is  no  longer  monotone  but  achieves  a 
maximum  value  at  a  finite  temperature.  This  again  puts  a  bound  on  \f„  but  now  the 
minimum  value  is  achieved  at  a  finite  temperature. 

This  predicts  a  maximum  possible  amount  of  blowing  from  behind,  that  is  a 
minimum  possible  value  of  A/„.  if  the  Lewis  number  is  less  than  one  and  the  reaction 
constant  it  is  less  than  or  equal  to  one  so  that  the  rate  of  reaction  approaches  zero  as 
the  temperature  becomes  infinite.  This  lower  bound  on  .V/„  is  surprisingly  independent 
of  the  composition  of  the  fresh  mixture.  The  maximum  possible  value  of  jA/ol  is 
attained  at  a  final  burned  gas  temperature  of  1/(1  -  u)  and  has  a  value  of 

IMiL.  =  Vf-C  -  «r"cxp(/i  -  1),.(1  -  L)  (41) 

This  means  that  as  the  blowing  from  behind  the  flame  is  increased,  extinction 
must  occur  at  or  before  the  value  of  the  blowing  given  by  (41)  if  both  L  and  n 
are  less  than  one. 

The  physical  reason  for  the  restriction  on  n  is  that  vvhen  11  is  less  than  one.  the 
reaction  rate  eventually  begins  to  fall  as  the  temperature  rises.  For  large  enough 
temperatures  the  drop  in  overall  concentration  with  rising  temperature  overwhelms 
the  increase  in  the  concentration  of  sufficiently  energetic  molecults  as  represented  by 
the  exponential  term.  This  behavior  can  appear  for  LeV^s  numbers  less  than  one  '*j 
because  the  blowing  from  behind  increases  the  temperature  of  the  mixture,  a  fact 
that  has  been  confirmed  through  numerical  integration  of  the  governing  equations. 
Sufficiently  large  amounts  of  blowing  can  conceivably  raise  the  temperature  to  the 
point  where  the  reaction  rate  begins  to  fall  due  to  a  sufficiently  large  drop  in  density. 

This  effect  is  the  reason  for  the  restriction  (41).  This  type  of  extinction  mechanism 
cannot  be  studied  through  the  use  of  AEA  since  it  can  only  occur  when  the  final  gas 
temperatures  arc  of  the  same  order  ol  magnitude  as  the  activation  temperature. 

Here  we  have  taken  the  mass  flux  of  the  remote  opposed  flow  to  be  finite  and  yet 
in  the  AEA/boundary  layer  solutions  the  magnitude  of  the  mass  flux  appears  to  grow 
without  bound.  There  is  no  inconsistency,  however.  In  the  AEA  solutions  the  velocity 
scale  is  the  adiabatic  flame  speed  whereas  here  wc  have  used  a  velocity  scale  that 
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m>mi;iii/CN  ific  ciK;iI'n.'icm  of  t‘ic  rcaciion  term.  The  AEA  soluiions.  where  the  iei;ipera- 
lure  o''the  remote  reaeianis  !s  small,  can  be  recovered  by  rcsiricune  T,  lo  be  small  and 
takine  M„  lo  be  laree  compared  to  the  adiabatic  burning  rate.  Thi.;  is  possible  even 
ihiiugh  the  critical  \  aluc  vif  M„  approaches  zero  as  T,  goes  to  zero  because  the  critical 
mass  flux  is  (){^  -(T, ))  whereas  the  adiabatic  burning  rate  is  0{^  -;(T,):T, )  and  so  the 
critical  mass  flux  is  much  larger. 

Arrhenius  kinetics  are  not  an  essential  feature  fo  this  extinction  mechanism.  If  we 
use  a  reaction  rate  temperature  dependence  of  the  form  /IT)  where  /(O)  =  0, 
f{T)  —  Oas  T  -»  r..  f(T)  >  0  for  7"  >  0,  and /(T) continuous  then  we  can  derive 
the  result  that 


A/,-;  <  Lj\T,)l{\  -  L)  (42) 

when  the  Lewis  number  L  is  less  than  one.  These  restrictions  on  f(T)  are  all  physically 
realistic  in  that  we  e.xpcct  essentially  zero  rate  of  reaction  for  small  temperatures, 
diminishing  rate  of  reaction  for  very  high  temperatures  since  the  fraction  of  sufTiacntly 
energetic  molecules  can  become  no  higher  than  one  and  the  concentration  ol  reactants 
must  decrease  wiih  temperature,  and  the  rate  of  reaction  must  always  be  positive. 
Since  the  physically  realistic  restrictions  imposed  on  /fT)  imply  that  f{T)  has  a 
ma.ximum  value,  we  can  therefore  conclude  that  there  is  a  maximum  amount  of 
blowing  that  a  flame  with  Lewis  number  less  than  one  can  withstand  and  this 
maximum  is  no  greater  than  that  predicted  by  (42). 

Although  an  extinction  mechanism  like  this  is  predicted  mathematically,  the 
question  remains  as  to  whether  or  not  it  is  physically  accessible.  In  addition,  it  is 
possible  that  the  failure  of  a  physically  realistic  flat-flame  similarity  solution  to  exist 
only  indicates  that  the  solution  bifurcates  or  jumps  to  a  nonplanar  solution.  A 
complete  answer  to  these  questions  would  require  the  study  of  a  complete  reaction 
mechanism,  large  scale  numerical  simulations,  developing  appropriate  existence  and 
uniqueness  theorems,  experimental  investigation,  or  some  combination.  Such  a  task 
is  beyond  the  scope  of  this  work. 


CONCLUDING  REMARKS 

Here  we  have  shown  that  there  exists  a  potential  mechanism  for  the  abrupt  extinction 
of  acrodynamically  strained  premixed  flames  with  Lewis  numbers  less  than  one.  The 
surprising  Icature  of  this  mechanism  is  that  instead  of  the  usual  notion  of  a  flame 
extinguishing  by  becoming  too  cold,  here  flames  can  extinguish  by  becoming  too  hot. 

The  physics  of  this  mechanism  is  surprisingly  simple.  When  the  Lewis  number  of  a 
combustible  mixture  is  less  than  one.  aerodynamic  straining  of  the  flame  will  cause  a  rise 
ill  flaivtc  tciripeiatuiC.  Eventually  the  rale  of  icaciion  must  fall  due  lo  the  fact  that  the 
concentration  of  active  species  is  inversely  proportional  to  the  temperature  (if  the 
pressure  is  fixed)  and  the  fraction  of  sufficiently  energetic  collisions  can  be  no  greater 
than  one.  This  falling  reaction  rate  with  increased  biow'ing  can  extingui.sh  (he  flame. 

The  question  remains  as  to  whether  or  not  this  extinction  mechanism  is  accessible 
to  physical  flames.  Perhaps  a  simple  experiment  can  reveal  the  answer. 
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ABSTRACT 

The  far  field  of  a  stretched  premixed  flame  with  A-B^C 
kinetics  is  studied.  Each  of  the  reaction  rates  are  assumed  to  be 
proportional  to  the  mass  fraction  of  the  reactants  and  have  a  tempera¬ 
ture  dependence  that  is  continuous,  positive,  and  vanishes  at  both 
infinite  and  zero  absolute  temperature.  Arrhenius  reaction  rate 
functions  are  examples  that  satisfy  these  requirements.  Bounds  are 
uncovered  for  tn-  ^mount  of  blowing  from  behind  the  flame  that  plane 
premixed  flames  can  withstand  for  all  Lewis  numbers  when  the  interme¬ 
diate  species  are  lighter  than  the  deficient  reactant.  When  the 
intermediate  products  are  heavier  than  the  deficient  reactant,  bounds 
on  the  blowing  are  found  when  the  Lewis  number  of  the  deficient 
reactant  is  less  than  one. 

INTRODUCTION 

Through  the  analysis  of  the  far  field  of  adiabatic  stretched 
premixed  flames  with  A-B  kinetics,  we  have  been  able  to  show  that 
premixed  flames  with  Lewis  numbers  less  than  one  will  be  extinguished 
given  sufficiently  strong  blowing  from  behind  [1].  This  is  inferred 
from  the  analysis  in  that  the  mass  fraction  of  combustible  becomes 
negative  in  the  far  field  if  the  opposed  flow  is  above  a  MXimum 
allowable  level.  The  activation  energy  asymptotic  (AEA)  analyses  of 
this  case  [2-7]  do  not  show  true  extinction.  What  we  mean  by  true 
extinction  is  the  failure  of  existence  for  a  flame  structure  for 
sufficiently  strong  opposed  flow. 

In  the  previous  work  [1]  it  was  unclear  as  to  what  effect 
dissociation  would  have  on  the  results.  Here  we  will  investigate  a 
kinetic  mechanism  that  incorporates  dissociation.  The  simplest 
possible  reaction  mechanism  that  would  include  dissociation  is  A«B. 
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It  seems  unJ  Ikely  that  the  analysis  of  such  a  mochariism  would  revea! 
anything  of  interest  in  that  the  final  equilibrium  state  would  be 
composed  of  major  species  only  without  any  intermediates.  A  more 
plausible  model  is  the  A> B< C  mechanism  where  a  represents  intermediate 
products  and  the  final  equilibrium  state  is  dominated  by  intermediate 
and  final  products. 

The  plan  of  this  paper  is  to  first  analyse  the  far  field  of 
the  stretched  A-B<C  flame  so  that  the  case  without  dissociation  is 
fully  documented.  Once  the  non-dissociating  model  is  fully  explored 
we  shall  include  the  effects  of  dissociation. 

MATHEMATICAL  MODEL 

The  equations  that  model  two-dimensional,  low  Mach  number 
combustion  waves  under  the  assumption  of  calorif ically  and  thermally 


perfect  gases  with  constant  transport  properties  are 

a (pu)/ax+a (pv)/ay»o  (i) 

a  (PUT) /3X+3  (P VT) /ay-72T+qj^Yf  (T)  +q2^^2  -<I2 ^ ^ 3  (T)  ( 2 ) 

3(puY)/3x+a(pvY)/ay-(LA)‘^^^y-yfj^(T)  (3) 

a  (puX)/ax+3  (pvX)/ay-{LB)  "-^^^X+Yf  ^  (T)  -Xf^  (T)  +Zf  ^  (T)  (4) 

2-1-X-Y  (5) 

puau/ax+pvau/ay— ap/ax+p(7^u-*'{a  (au/ax+3v/ay)/ax)/3)  (6) 

puav/ax+pvav/ay— ap/ay+P(v^v+(a  (au/ax-*-av/ay)/ay)/3)  (?) 

O-l/T.  (8) 


u  and  V  are  the  non-dimensionalized  x  and  y  components  of  fluid 
velocity,  T  is  the  non-dimensional  temperature,  p  is  the  non- 
dimesional  density,  p  is  the  non-dimensional  pressure,  P  is  the 
Prandtl  number  and  X,  Y  and  Z  are  the  mass  fractions  of  the  intermed¬ 
iate  product,  the  reactant  and  the  final  product  respectively.  These 
equations  are  continuity  (1),  energy  balance  (2),  species  balance  (3- 
?) ,  X  and  y  momentum  balance  (6-7)  and  equation  of  state  (8). 

Proceeding  as  in  [1],  this  system  can  be  simplified  for 
plana  flaaMS  in  a  nun-uniform  flow  field.  The  resulting  equations 
under  the  limit  of  large  heat  release  are 


KdT/dx-d^T/dx^+q^Yf - (T) (T) -qj ( 1-X-Y) f 3 (T) 

(9) 

MdY/dx- ( l/L*) d^Y/dx*-Yf 3 (T) 

\ 

(10) 

MdX/dx- ( I/Lb) d^X/dx*+Yf 3 (T) -Xf ^ (T) + ( 1-X-Y) f 3 (T) 

(11) 

Pd^  (TdM/dx)/dx^-Md(TdM/dx)/dx+T(dM/«ix)  ^-0. 

(12) 

H  is  the  mass  flux,  pu.  without  the  limit  of  large  heat  release 
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^^'na.W"^reinoT:e  reactants  side  of  equation  (12) 

could  be  a  non-zero  constant. 

Appropriate  boundary  conditions  are 

X(-'-)»f3(Tf)/(f2(Tf)+f3(Tf))  ,  Y(+-)-0  (13) 

where  Tf  is  the  temperature  of  the  remote  products,  a  quantity  that 
must  be  found  as  part  of  the  analysis.  In  addition,  a  boundary 
condition  is  needed  to  force  the  balance  in  the  remote  reactants  to  be 
one  between  conduction  and  convection  without  appreciable  reaction. 
Such  a  condition  is 

(d2T/dx2)/(M  dT/dx)  -  1  as  x  -  — .  (14) 

THE  REMOTE  BURNED  STATE 

Since  the  system  (9-13)  is  autonomous,  we  can  transform  the 
problem  into  a  phase  space  problem  on  a  finite  domain  with  independent 


variable  T  through  the  definition 

z  -  dT/dx.  ^  (15) 

The  resulting  system  of  equations  is 

M  -  dz/dT  (qiYfi(T)  +  q2Xf2(T)  -  q2  ( 1-X-Y)  f  3  (T) ) /z  (16) 

MdY/dT  -  (1/La)  d(2dY/dT)/dT  -  Yfi(T)/z  (17) 

MdX/dT  -  (1/Lb)  d(2dX/dT)/dT  + 

(Yfi(T)  -  Xf2{T)  +  (1-X-Y) f3(T))/2  (18) 

?zd(zd(T2dM/dT)/dT)/dT  -  M2d(TzdM/dT)/dT  +  Tz2{dM/dT)2  .  0  (19) 

and  the  boundary  conditions  at  X"-*-»  become 
Z(Tf)-0,  X(Tf)-f3(Tf)/(£2{Tf)+f3(Tf))-Xf, 

Y(Tf)-0,  M(Tf)-Mo.  (20) 

Near  T*Tf  we  can  write 
z-<r,  Y-<*,  x-Xf-*-<f 

T  -  Tf  M  -  Mo  Sn.  (21) 

The  resulting  system  is 

Mo  -  df/dr  +  qi^fi(Tf)  + 

q2C^A(Tf)  +  *f3(Tf)  +  {(£2(Tf)+f3(Tf))]  (22) 

Mo  :  d#/dr  -  (1/La)  r  cl(r  d^/dr)/dr  -  ^fi(Tf)  (23) 

Mo  r  d£/df  -  (1/Lb)  f  <i(f  d«/dr)/dr  ^fi(Tf)  - 

f  (f2(Tf)-^f3{Tf))  -  ^f3{Tf)  -  rA(Tf)  (24) 

and  another  equation  for  m  that  is  of  no  consequence  in  the  remainder 
of  the  analysis.  The  function  a  is  defined  by 

A(T)  -  (£3  d£2/dT  -  f2  df3/dT)/(f2  ♦  £3)-  (25) 

The  boundary  conditions  are 

r(0)-^(0)-? (0)-0. 


(26) 
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The  singular  point  described  by  the  system  (22-24,  26)  is 
inherently  non-linear,  but  fortunately  it  admits  solutions  of  the  form 
0=ar.y=br,f»Cr.  (27) 

Substitution  into  (22-24)  results  in 

c  =  [LaMo  -  .'{LaMo)^  t  4LACi(Tf)]/2  (23) 

and  the  set  of  algebraic  equations 

a(qifi(Tf)+q2f3(Tf))  -*•  bq2  (  f  2  (Tf )  ■►f  3  (Tf )  ) 

=.  MqC  -  c2  -  qjA  (Tf)  (29) 

a(fi(Tf)-f3(Tf))  +  b(c2/LB  -  ”0=  “  f2(Tf)  “  f 3 (Tf ) )  -  A (Tf )  (30) 

that  are  trivially  solved  for  a  and  b  so  that  a,b  and  c  are  now  known 
as  functions  of  Mq  and  Tf.  There  is  another  solution  for  c  with  the 
positive  root  in  (28).  but  that  solution  is  not  consistent  with  the 
fact  that  T  approaches  a  constant  as  x  goes  to  infinity. 

THE  REMOTE  SINGULAR  POINT  FOR  NO  DISSOCIATION 

The  solution  fails  to  be  physically  realistic  when  a  is 
positive  since  that  would  correspond  to  a  negative  mass  fraction  for  'i 
near  the  singular  point.  The  case  where  a*0  therefore  defines  a 
critical  condition. 

The  specification  of  the  critical  condition  is  very  simple 
in  the  case  of  no  dissociation.  When  dissociation  is  neglected  we 
must  take  f3(T)-0  which  implies  that  a-0.  Setting  a  equal  to  zero  and 
eliminating  b  from  (2'Ks30)  gives  rise  to 

(MoC(l-LA)-fi(Tf)LA) (MoC(La/Lb  -  1)+La/Lb  f l (Tf ) -f2 (Tf ) ) -0.  (31) 

Substituting  (28)  for  c  gives  the  two  critical  conditions 

Mo^c  “  LAfi(Tf)/(l-La)  (32) 

and 

Mo^c  *  (^2(Tf)-LA/LB  f l (Tf ) ) V [ La ( f2 (Tf ) -f 1 (Tf ) ) (La/LB-1) } .  (33) 

Since  the  Lewis  numbers  and  reaction  rates  are  always  positive, 
condition  (32)  can  only  be  attained  when  La  <  1.  In  addition,  if  we 
only  consider  the  case  where  f2(T)  >  fi(T)  throughout  the  range  of 
possible  final  product  temperatures  Tf  so  that  product  formation  is  a 
faster  reaction  than  the  breakdown  of  the  reactants  (usually  a  good 
assumption  wince  the  breakdown  of  the  relatively  stable  reactants  is 
governed  by  a  fairly  high  activation  temperature  mechanism  and  the 
formation  of  product  from  intermediates  is  typically  through  reactions 
with  small  activation  temperatures) ,  then  condition  (33)  can  only 
apply  if  La  >  Lb-  This  is  a  very  common  occurence  since  it  requires, 
roughly  speaking,  the  mean  molecular  weight  of  the  intermediate 
species  to  be  smaller  than  the  molecular  weight  of  the  deficient 


71 


reactant.  The  only  potential  mechanism  that  we  can  think  of  that  does 
not  conform  to  this  condition  is  lean  hydrogen  oxidation  where  the 
deficient  reactant  has  molecular  weight  2  and  the  intermediate  species 
are  H,  OH  and  O.  It  is  likely  that  the  mean  molecular  weight  of  the 
intermediates  will  be  greater  than  2  in  this  case. 

Condition  (32)  is  Identical  to  the  critical  condition  for 
the  A->B  mechanism  [1]. 

An  example  where  <  Lx  and  Lx  >  1  is  given  in  Figure  1. 
This  situation  occurs  when  the  intermediate  products  are  ligther  than 
the  deficient  reactant  and  the  deficient  reactant  is  heavier  than  the 
bulk  gas.  In  this  case  only  .(33)  defines  a  limiting  condition.  In 
this  and  all  subsequent  figures,  the  shaded  regions  are  regions  in  the 
parameter  plane  where  the  remote  reactant  mass  flux  is  mathematically 
negative  and  hence  is  of  no  physical  significance. 
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Fig.  1:  Critical  remote  mass  flux  versus  flame  temperature. 

Lx-1.2,  Lb-. 8,  Qi— 1,  q2“2,  f i-exp(-l/T)/T,  f 2-exp (-. 1/T)/T 

Figure  2  shows  a  case  when  Lx  <  1  end  Lx  <  Lg*  Since  this 
case  only  applies  to  the  situation  where  the  deficient  reactant  is 
both  lighter  than  the  intermediate  products  and  the  bulk  gas,  it  seems 
likely  that  it  may  only  apply  to  lean  hydrogen/oxygen  or  lean  hydro¬ 
gen/oxygen/diluent  mixtures.  In  this  case  only  the  condition  (32) 
applies. 

Figure  3  shows  a  case  where  Lb  <  ^x  <  1.  Such  a  situation 
is  likely  for  lean  methane/air  flames,  for  example,  where  the  prin-* 
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clple  interinod latoa  aru  1  i'/hCer  than  the  dotlcionc  reactant  which  in 
turn  is  lighter  than  the  bull:  fluid.  In  this  case  critical  condi- 
tionn  are  given  by  both  (32)  and  (3J). 

In  tho  first  two  cases  it  is  clear  that  there  are  no  physically 
realizable  solutions  if  the  remote  mass  flux  in  the  burned  gas 
becomes  too  high.  In  the  last  case  shown  in  Figure  3,  the  mass 
fraction  of  the  remote  reactants  is  positive  if  conditions  are  such 
that  the 


Fig.  2:  Critical  remote  mass  flux  versus  flame  temperature. 

La-. 5,  Lb-. 7,  qi— 1,  q2-2.  f i-exp (-1/T)/T.  f2-exp(-. 1/T)/T 

flame  is  above  both  curves.  As  a  practical  matter,  however,  such 
solutions  are  probably  not  accessible  since  the  starting  conditions 
in  any  real  system  would  be  in  the  lover  domain  and  there  is  no  way 
to  enter  the  upper  domain  in  any  continuous  manner  except  through  the 
single  point  where  the  two  domains  meet.  The  "solution"  at  the  point 
where  the  two  solutions  meet  is  found  as  limit  of  problems  where  the 
remote  upstream  reactant  mass  fraction  goes  to  zero  and  hence  is  not 
truly  a  flame.  In  fact,  the  solid  curves  are  not  part  of  the  region 
of  physically  realistic  flame  solutions  because  they  correspond  to 
situations  where  the  remote  upstream  mass  fraction  of  reactant  is 
zero  and  hence  there  is  no  flame.  This  fact  is  found  chrough 
numerical  integration  of  the  governing  equations. 


t.\ 

THE  REMOTE  SIMCUIJkR  POINT  WITH  DISSOCIATION 

Whan  dlHsociatlon  la  ratainad,  thw  critical  Hq  is  fcurxl  by 
solving  tha  algabraic  aquatlor 

(CMo(  1-I^a)  -L^f  1  (Tf )  -q2A  (Tf ) )  (CMo(-1*Lx/Lb)  ♦f  1  (Tf )  L*/Lb 
-f2tT.)-r3(T/))  -  q2A(rf)(f2(Tf)^f3(Tf))(34) 
whare  C,  a  function  of  Mq.  !•  qivan  by  (28).  This  is  dons  aost  easily 
by  numarical  root  finding. 

First  ws  will  considar  tha  case  whara  Lg  <  and  >  1, 
such  as  a  laan  heavy  hydrocarbon/a I r  praaixad  flana,  with  weak 
dissociation. 


Figure  3:  Critical  ranote  mass  flux  versus  flame  temperature. 

La-. 8,  Lb-. 4,  qi— 1,  q2-2,  fi-exp(-l/T)/T,  f2-exp(-. 1/T)/T 

A  typical  parameter  plot  is  shown  in  Figure  4.  Comparing  this  plot  to 
the  corresponding  case  without  dissociation  (Figure  1)  we  see  that 
dissociation  decreases  the  maximum  amount  of  blowing  that  the  flame 
can  withstand. 

tn  the  case  where  La  <  Lb  <  1,  «•  shown  in  Figure  5,  the  maximum 
possible  amount  of,  blowing  is  increased  through  the  action  of  dissoc¬ 
iation.  As  mentioned  earlier,  such  an  ordering  of  the  Levis  numbers 
is  very  infrequent  in  applications.  One  of  the  only  possible  physical 
systems  that  night  have  such  an  ordering  is  lean  hydrogen/air  premixed 
flames. 

The  case  where  Lb  <  La  <  1  i*  shown  in  Figure  8.  This  manner  of 


n 


ramotc  burnad  A  ftp  op«nfl  up  through  which  colutlon  curvoo 

may  pass  and  tha  atrongor  tha  rata  of  dissociation,  tha  widar  tha  gap. 
It  ia  probabla, though  wa  hastan  to  add  that  wo  hava  not  chaekad  this 
raault  fully,  that  situation  can  ariaa  whom  naar  atoiehioaatric 
mlxturaa  will  axtlngulsh  with  suffieiantly  high  blowing  but  rich 
mlxturaa  will  net  with  so«m  pivotal  value  of  tha  ranota  aaoa  fraction 
of  raactant,  danotad  by  Ye,  that  saparatas  tha  two  caaas.  This 
phanomana  ia  shown  schsaatically  in  Figura  7.  As  tha  rata  of  dissoc' 
iation  is  incraasad,  tha  gap  should  widan  with  a  subsaquant  incraasa 
in  tha  pivotal  valua  of  tha  ranota  raactant  nass  flux.  With  a 
suffieiantly  high  rata  of  dissociation,  it  saaas  likaly  that  tha 
pivotal  valua  of  tha  nass  fraction  of  raactants  should  go  to  1  with 
furthar  incraasas  in  tha  rata  of  dissociation  rasulting  in  no  aMxinim 
rata  of  blowing. 


Figura  6:  Critical  nass  flux  versus  flana  tanparatura. 
La". 8,  Lb". 4,  qi"-i»  q2*2,  fi"axr(-i/T)/T 

f2-axp(-.l/T)/T,  fj-.OOl  axp(-2/T)/T 


CONCLUDING  REMARKS 

Kara  wa  hava  analytad  tha  response  of  a  plana  praaixad  flana  with 
saquantial  kinatics  to  a  non-unifom  flow  fiald.  Two  aaparata  casas 
wara  invast igatad;  ona  vhara  tha  sacond  stap  was  irravarsibla  and  ona 
whara  it  was  ravarsibla.  Tha  inclusion  of  ravarsibilty  in  tha  sacond 
stap  did  not  altar  tha  qualitative  flana  raaponaa  axoapt  in  tha  casa 
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